Order Parameter "Holes" and Theory of Microwave Conductivity in YBa 2 Cu 3 



7-5 



OS 

On 



3 



o 

CO 



C 

O 

o 



> 
o 
in 

o 

On 

ON 

-i— > 
o3 



i 

C 

O 

o 



Matthias H. Heftier"'* and Peter J. Hirschfeld 6 

a j4n?onne National Laboratory, 9100 South Cass Ave, Argonne, IL 
b Department of Physics, University of Florida, Gainesville, FL 32611, USA. 

We propose that the low temperature discrepancy between simple d-wave models of the microwave 
conductivity and existing experiments on single crystals of YBa2Cu307_« can be resolved by includ- 
ing the scattering of quasiparticles from "holes" of the order parameter at impurity sites. Within a 
framework proposed previously, we find in particular excellent agreement with data of Hosseini et 
al. on slightly overdoped YBa2Cu307_a samples over the entire temperature range down to about 
2-3K, and for a wide range of frequencies. Remaining discrepancies in the "universal" regime at 
very low temperatures are discussed. 



PACS Numbers: 74.25.Dh, 74.25.Nf, 74.62.Dh 

Introduction. One of the few areas of the high- 
temperature superconductivity problem where a well- 
defined theoretical description is available is the super- 
conducting state of the optimally doped materials. Al- 
though the normal state of the cuprates cannot be de- 
scribed in terms of weakly interacting Landau quasipar- 
ticles as in a simple metal, it appears that such states 
are well-defined below the critical temperature, Q due in 
part to the well-known collapse of the quasiparticle scat- 
tering rate. |||| This justifies the application of the BCS 
theory, with a dominant pairing component which is now 
widely agreed to correspond to d x 2_ y 2 symmetry. The 
successes of the (i-wave theory are many, and have been 
amply reviewed, j^] The agreement of experiment and 
theory is particularly impressive in the YP^CuaOy-a 
(YBCO) system, where steady advances in the quality 
of single crystal growth in the past several years have 
allowed the clean d-wave superconducting state to be 
probed. 

One qualitative discrepancy that remains in the com- 
parison of d-wave theory with bulk measurements prob- 
ing the ab plane concerns the low-temperature microwave 
conductivity cr(f2, T). Measurements on yttria-stabilized 
zirconia crucible-grown, nominally pure, YBCO crystals 
at low frequencies f2 of order a few GHz showed an appar- 
ently linear temperature dependence, a ~ T for tempera- 
tures T sufficiently far below a large peak at 30-40K. 
Theories of the dynamical conductivity in the d-wave 
state have been successful in explaining roughly the 
size and position, as well as the frequency and disor- 
der dependence of the peak. However, at low tempera- 
tures and small microwave frequencies, where transport 
is dominated by elastic scattering and the theory should 
be simplest, the theoretical prediction is a ~ T 2 , not T. 

A pure linear temperature dependence turns out to be 
very difficult to obtain within the framework of the BCS 
theory, however. In the original work of Bonn et al., Q it 
was suggested that er ~ T is a natural result at low tem- 
peratures since the assumed Drude form of the conduc- 
tivity a ~ (e 2 /m)n qp (T)T(T) yields a ~ T if the quasi- 
particle density varies as n qp ~ T as expected for a ci- 
wave superconductor at low energies, and if the effective 
quasiparticle scattering time r(T) saturates at low tem- 



peratures, as in a normal metal. While the assumption of 
a Drude form of the conductivity was supported by mi- 
croscopic analysis, it was shown that pair correlations 
in the usual impurity scattering models generically lead 
to strong temperature dependence of the scattering time, 
e.g. r(T) - T (unitarity limit) or r(T) ~ 1/T (Born 
limit). Attempts to resolve this problem by choosing 
intermediate scattering strengths || have not provided 
obviously better results. 

Ultimately, the problem can be traced back to the fact 
that the BCS Green's function has quasiparticle poles 
at uj = ±I?k, where = v^k + A 2 ,, £k is the single- 
particle band energy, and At = Ao cos 2cf> is the d-wave 
order parameter (OP) over an isotropic 2D Fermi sur- 
face. This leads naturally to an analytic expansion of the 

— > conductivity in powers of T 2 , regardless of band 
structure or phase shift (though the range where this 
expansion is possible can vary dramatically). The lead- 
ing term in the d-wave case is the "universal" (disorder- 
independent) constant er(f2 — * 0,T — ► 0) = o"oo = 
ne 2 /(m7rA ), || hence the T 2 result obtained for the 
subleading term in Ref. Q. While initially microwave 
cavity experiments did not have the resolution required 
to determine this constant, its analog in the case of ther- 
mal conductivity was confirmed recently by Taillefer et 
al. @ Whether the same discrepancy regarding the sub- 
leading temperature corrections exists in the case of ther- 
mal currents is not known, since the electronic corrections 
are obscured by the phonon contribution. 

The microwave conductivity impasse was broken very 
recently with data on extremely pure BaZr03 crucible- 
grown YBa2Cu3 06.99 single crystals by Hosseini et al. 
flllf In addition to observing the increase in peak height 
and decrease in peak temperature predicted by the Drude 
picture for purer samples, it was noted that the very high- 
resolution data did not support the a(fl — > 0, T) ~ T re- 
sult anticipated on the basis of the earlier measurements. 
Rather, the T dependence of the lowest frequency (~ 1 
GHz) data was slightly sublinear, crossing over to super- 
linear dependence above a crossover frequency of about 
13 GHz. This unusual behavior has led us to reconsider 
a proposal we advanced recently, jl2| that the standard 
d-wave theory of impurity scattering should be modified, 
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particularly in the case of transport properties, to include 
the effects of order parameter "holes" near the impurity 
sites. [[l2|-[l4[ This leads to additional scattering at low 
energies due to the formation of an additional scatter- 
ing resonance. |[13] A recent calculation on the spatial 
structure of the local density of states around impuri- 
ties sites has used similar methods. [ [j"5| In this paper we 
argue that, with the proper inclusion of the corrections 
due to OP suppressions, nearly all features of the low-T 
microwave conductivity data can be understood quanti- 
tatively. The exception is the value of the residual O — > 
conductivity, which we discuss in some detail. 

Impurity t-matrix. We assume in our approach jl^] 
that the bare impurity can be described by a ^-function 
scattering potential, [/(R — Ri mp ) = [/o<5(R — Rimp)T3, 
where the Ti are the Pauli matrices in particle-hole space. 
We then argue further |fl2f that for bulk transport prop- 
erties the only essential features of the order parameter 
suppression induced by this potential are that it retains 
the roughly the same symmetry as the background or- 
der parameter (induced subdominant pair components 
can be included but we deem them less important) and 
deviates substantially from the bulk value over a range 
of order the coherence length £o or smaller around the 
impurity site. The smallness of £o in the cuprates for 
T <C T c is then used to justify our replacement of the 
order parameter fluctuation by a point like potential, 
<5Ak(q = 0) <5(R — Ri mp ), where Ri mp is the impurity 
site and we take (5 Ak(q = 0) = Sd cos 2<\>. In Ref. [[l2| the 
amplitude Sd is then determined self consistently through 
the BCS gap equation for a single impurity. Both the gap 
equation and the conductivity, calculated now by stan- 
dard impurity-averaging methods, then contain the effec- 
tive impurity potential U(p, p') = (UqTt, + 5dCOs2<j)Tij . 

Thus we explicitly account for the fact that electrons 
moving near the impurity feel an effective one-body po- 
tential due not only to the bare impurity but to the order 
parameter fluctuation around it. 

The impurity t-matrix is given as usual by (in matrix 
notation) T = U + UGqT, with U defined above. In 
the usual "dirty d-wave" theory, the t-matrix is taken 
independent of momentum, T = T(oj), for the case of 
isotropic scatterers. It becomes momentum-dependent 
in the current theory with d-wave OP suppression. 

The solution for the t-matrix at q = p' — p = in the 
present ansatz for a single impurity may be written 



eraged G via the Dyson equation, G 
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where go and g 2 are the components of the momentum 
integrated Green function, g = (1/2) ^ fc TrG(k, w) and 
g 2 = (1/2) J2k ^ r Tl cos20G(k, uj). Again, subleading 
OP contributions have been neglected. |1J] The disorder- 
averaged self-energy is now defined in the limit of a 
density n.j of independent impurities to be S(k, uj) = 
rtiTk(w), and determined self-consistently with the av- 



AfcTi — S(k, u>) = u) — S,kT~3 — AfcTi. The first term in Eq. 
[j] is the usual dirty d-wave theory result for arbitrary 
scattering phase shift So = — cot _1 (l/(iVo7r[/o)); yielding 
in particular a resonance at u> = in the unitarity limit 
Uo — * oo. 

The denominator in the second term, due to OP scat- 
tering, leads to a similar resonance, at a position to be 
numerically determined. Introducing Cf = (irSdNo)^ 1 we 
estimated in Ref. the quantity uj/ Aq — [— 2/ir — c/] = 
—Cf which determines the position of the off-diagonal 
resonance to be c/ ~ —0.16 as Uo — > oo. It is impor- 
tant to note that this determination was based on a solu- 
tion to the one-impurity problem neglecting the average 
suppression of the gap due to disorder, and without ac- 
counting for induced subdominant pair components. We 
expect that \Sd\ in a complete theory to be somewhat 
smaller, and for the moment take c/ to be a free param- 
eter. The main point, however, is that some additional 
spectral weight is shifted from the gap edge down to low 
and intermediate frequencies ]T^ ] when compared to the 
usual d-wave model. We now investigate how this affects 
the microwave conductivity and compare to the experi- 
mental data of Hosseini et al. |ll[] 

Microwave conductivity. In Ref. |Q, we discussed the 
51 — ► conductivity and compared with the expected 
a ~ T behavior; at that time it was noted that no pure T- 
linear behavior could be identified in the theory, although 
the deviations from the theory without off-diagonal scat- 
tering were pronounced and of the right qualitative size 
and sign. The data of Hosseini et al. |ll]] now make it 
clear that a subtle crossover in the temperature depen- 
dence is taking place in the l-20GHz regime. We there- 
fore calculate the full frequency-dependent conductivity 
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where £± = ± y Cj ± — A £ ± , the subscripts ± indicate 

evaluation at w±iO + , and primed quantities are evaluated 
at u> — O. The first term in Eq. ^ gives rise to the "uni- 
versal" T — > conductivity croo, while the second term 
determines the leading temperature corrections. As in 
Ref. , quasiparticle states are broadened in all calcula- 
tions by an inelastic scattering term 1 /ri„ e ; (T) calculated 
from a model of spin-fluctuation scattering as described 
in Ref. @. 

In Fig. |l| we show the best fit obtained to the low- 
frequency data of Ref. jllj , with normal state scattering 
rate T = 0.0005T C in the unitarity scattering limit, and 
take Cf ~ —.95 and A /T c = 3. We note that the the- 
ory indeed reproduces the slightly sublinear behavior at 
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low frequencies, and crosses over at around 13GHz to 
a super linear behavior and eventually to T 2 . At low T 
the theory curves are consistently below the experimen- 
tal data, but the relative discrepancy is not large, with 
the exception of the very low temperature regime below 
2-3K. 
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FIG. 1. Normalized microwave conductivity for the fre- 
quencies of Ref. jff[. Solid lines: theory with T/T c = 0.0005, 
c f = -.95, A /Tc = 3. Dashed line: same T, / = 2.25 GHz, 
without OP scattering. Dot-dashed line: T = .0014, / = 2.25 
GHz, without OP scattering. 

In Fig. [2] we show the frequency dependence of the 
theory for the same parameters at experimental temper- 
atures; here the agreement is excellent. As noted by Hos- 
seini et al., the most remarkable aspect of the data is the 
extremely weak temperature dependence of the width in 
frequency of the residual conductivity peak at low tem- 
peratures; this is well reproduced by the theory. In the 
insert to Fig. |^, we plot the experimental determination 
of this width, called l/r(T), obtained by a fit to a Drudc 
(Lorentzian) spectrum. (jl]J However, in the microscopic 
model we propose, the relaxation rate 1/t(lu) for nodal 
quasiparticles is frequency-dependent. Furthermore, we 
find that the momentum dependence of the off-diagonal 
impurity self-energy prevents a derivation similar to Ref. 
providing a microscopic justification for a Drude or 
"two-fluid" analysis. Thus I/r(T) has no well-defined 
microscopic meaning, but for comparison's sake we have 
also performed a Drude fit to the theoretical (j(uj, T), the 
result of which is shown in the insert to Fig. 0. We 
note that the deep minimum in l/r(T) found without 
the OPS |?| has been nearly eliminated, but the scat- 
tering rate still rises at the lowest temperatures due to 
the increase of the elastic scattering rate in the unitar- 
ity limit (see inset). However, this rise is somewhat 
spurious as the conductivity spectrum becomes distinctly 



non-Lorentzian below 7 ~ \/TAo ~ 2-3K, indicating a 
further breakdown of the Drude interpretation. We note 
that in the present theory the minimum of l/r(T) occurs 
at an energy of roughly l/r mi „ = (r 3 A ) 1/4 . 
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FIG. 2. Normalized a(Q) at various T. The open sym- 
bols are experimental data from Ref. Jn|, the solid sym- 
bols are results from our theory, with F/T c = 0.0005 and 
Cf — —0.95. Insert: Effective scattering rate t(T c )/t(T) vs. 
T. The solid line corresponds to theory with parameters as 
above, the dashed line is the theory without OP scattering. 
The symbols are data from Ref. |fl|. 

Finally, in Fig. |3] we show the behavior of the con- 
ductivity over the entire temperature range, where the 
inelastic scattering plays a much more crucial role. It is 
remarkable that with the same parameters used to de- 
scribe the entire frequency range at low temperatures, 
the peak position is also well described. However, the 
theoretical curves rise faster than the experimental data 
below T c and fall faster than experiment below the peak 
temperature. At first glance this seems to be a shortcom- 
ing of our treatment of the inelastic scattering following 
Ref. On the other hand, the inset of Fig. || shows 

that the scattering rate extracted from the experimen- 
tal a-axis penetration depth and conductivity data via 
1/t(T) ~ n"(T)/q-(T) is in excellent agreement with the 
result of Ref. lfll in the T range dominated by inelastic 
scattering (Here, n"(T) = I - \ 2 a (T = 0)/A|(T) is the 
normal fluid fraction). We therefore believe that the ne- 
glect of band structure effects in the calculation of the 
density of states is the major source of the discrepancy 
at higher temperatures. 

Conclusions. We have shown that the inclusion of the 
order parameter suppression <5Ak at impurity sites dra- 
matically improves the fit of the weak-coupling BCS d- 
wave theory of the microwave conductivity to recent ex- 
periments on very pure samples of YBa2Cu307_a . The 
excellent fit over a wide range of temperatures and fre- 
quencies is found for an impurity scattering rate of ap- 
proximately r/T c = 0.0005, confirming that the sam- 
ples measured by Hosseini et al. are significantly purer 
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than the previous generation of yttria-stabilized zirconia 
YBa 2 Cu306.95 crystals. Although our treatment of SA^ 
is crude, it captures the basic physics of the curate su- 
perconductors, namely that the bulk order parameter is 
significantly suppressed over a very short length scale. 
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FIG. 3. o-(T) at various external frequency £1. The sym- 
bols are experimental data (Ref. pJ| , with the sharp peak 
at T c truncated). The solid lines are results of our theory 
with F /T c = .0005, Cf = —.95. Inset: Comparison of inelastic 
scattering rate after Ref. JhJ and the experimental ratio of 
normal fluid fraction and conductivity at 1.14 GHz. For the 
considered temperature range the impurity contributions to 
the scattering rate are unimportant. 

There are two sources of concern which suggest that 
the current approach needs to be refined still further. 
The first is that we were forced to treat the position Cf 
of the off-diagonal resonance as a fit parameter rather 
than determining it self-consistently, as originally pro- 
posed in jl2). We have suggested that the facts that i) 
the value determined from the fit is quite close to the 
self-consistently determined value (c/ ~ —0.8 with weak 
T dependence in the strong scattering limit), and ii) the 
approximations made, namely neglecting the effect of dis- 
order on the average gap value in the determination of c/ 
as well as induced subdominant pair components, tend 
in the right direction to account for the discrepancy, but 
we have no estimate of the magnitude of this error. Both 
a more complete, self-consistent calculation within our 
framework, as well as a full calculation of the one-electron 
spectral function in the framework of the Bogoliubov-de 
Gennes equation are desirable. 

The second problematic aspect is the limiting value 

and behavior of the conductivity for T<2 — 3K. As seen 
in Figure 1, the T — > value appears to be 3-5 times 
the calculated a m /a{T c ) = r" 1 (T c )/(7rA ). We have 
considered several kinds of corrections to the theory as 
described herein which tend to increase coo, or the ex- 
trapolated a(T — > 0). The first are nonlocal corrections, 
which are negligible in the geometry considered by Hos- 
seini et al. The second are deviations from the unitarity 



limit, which have been discussed by in Ref. ||. We find, 
as these authors did, that such deviations rapidly lead to 
increasing positive curvature in a vs. T and are never 
consistent with the measured cr(T) when a(T — > 0) is 
large enough to resemble experiment. More likely sources 
of error include the neglect of vertex corrections due to 
anisotropic impurity scattering or Fermi liquid correc- 
tions, which may change the residual conductivity by a 
factor of order unity. |18j The final possibility is that a 
few twin boundaries are present within the skin depth of 
the experiment. Twins have been shown to dramatically 
increase the residual conductivity in poorer samples, 
and lead to an "extrinsic" contribution to the transport 
which is poorly understood at present. 
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